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Abstract 

In a recent paper [I. Baldea and H. Koppel, Phys. Rev. B 78, 115315 (2008)], we showed that 
a variational approach [P. Delaney and J. C. Greer, Phys. Rev. Lett. 93, 036805 (2004)] proposed 
to compute the electron transport through molecules, which is based on boundary constraints of 
the Wigner function, is unable to correctly describe the zero-bias conductance of the simplest 
uncorrelated and correlated systems. In the present paper, we extend our previous analysis of 
the linear response limit of that approach, by considering, instead of the Wigner function, general 
constraints. We demonstrate that, if, as usual in transport theories, the quasi-particle distributions 
in electrodes are constrained, this method yields the completely unphysical result that the zero-bias 
conductance vanishes. Therefore, we conclude that the variational approach itself is defective. 

PACS numbers: 73.63.-b, 73.23.-b 
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I. INTRODUCTION 



In spite of significant advances, the field of molecular electronics continues to remain con- 
fronted with problems unsatisfactorily clarified both from the experimental and the theoret- 
ical sides. One major issue is the large discrepancy, often of several orders of magnitude,- - - 
between the measured currents through molecules and values calculated theoretically. Open 
questions in electron transport through nanostructures, which remain unresolved after many 
years of research, appear also in other related areas, e. g., transport through semi conducing 
quantum dots 4 ' 5 

The simplest theoretical calculations of the transport in nanoscopic/molecular systems 
rely upon the Landauer formalism,-^ which ignores electron correlations. The latter are 
accounted for in more elaborated approaches, like those based on nonequilibrium Green 
functions (NEGFs), 8 ' 9 time- dependent density matrix renormalization group (DMRG),— ~— 
and numerical renormalization group (NRG).—. In the most popular ab initio calculations 
in molecular electronics, the NEGF method is combined with calculations based on the 
density functional theory (DFT).— ~— Other approaches are based on various many-body 
schemes.— _ — Although the role of uncontrollable experimental factors should not be under- 
estimated (compare, for example, the contradictory results for one and the same system^"—), 
the ability to predict currents comparable with measured values continues to represent an 
important challenge for theory. 
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seemed to 



Out of the existing theoretical methods, that proposed and used in Refs. 
be very appealing, because it claimed to correctly reproduce currents experimentally mea- 
sured through correlated molecules and yielded a few other plausible results. It represents 
a variational approach with constraints imposed on the Wigner function at the device- 
electrode interfaces. The latter aspect is original: in usual transport theories, the electron 
distributions are constrained to be the Fermi functions. In a recent paper,— we inquired the 
validity of this method within the realm of theory and demonstrated that, for the simplest 
uncorrelated and correlated systems, its results are completely unphysical. However, because 
our critique primarily envisaged the imposed boundary conditions, one might still hope that 
this approach can be remedied by imposing other constraints. To investigate whether this 
is the case or not is the main purpose of our present work. 

The remaining part of the paper is organized in the following manner. The variational 
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24, will be de- 



approach considered here, which generalizes the approach proposed in Ref. 
scribed in Sect. [TTJ Its linear response limit will be considered in the next Sect. IIHI As 
an important point of the present work, which is emphasized in Sect. IIVI we argue that, 



contrary to what was claimed in Ref. |24j, strong electron correlations in the device do not 
preclude to impose usual boundary conditions, namely, to constrain the single-electron dis- 
tributions in electrodes to be the Fermi distributions. The main result deduced in Sect. 
HVlis that, even with these usual boundary conditions, the variational approach is inappro- 
priate: it unphysically predicts a zero-bias conductance that vanishes in situations where 
this is definitely not the case (e. g., uncorrelated systems, Coulomb blockade peaks, Kondo 
plateau). Conclusions are presented in Sect. [V] 



II. DESCRIPTION OF THE VARIATIONAL APPROACH 

Some years ago, Delaney and Greer (DG) proposed a method to calculate the steady-state 
electric transport through correlated molecules connected to two electrodes (two-terminal 
transport set-up) within a many-body approach.— It relies upon many-body calculations for 
a finite cluster consisting of the (nano)device (the molecule of interest) and (small) parts of 
electrodes with and without voltage bias. Their approach can be summarized as follows. A 
constant voltage bias, which is described by a term W in the Hamiltonian, drives the cluster 
at zero temperature (T = 0) from its ground state \l/ (H\^ ) = £ |\I> )) to a steady state 
\l/ characterized by a constant current flow. 

(i) The steady state wave function \l/ is determined by minimizing the total energy 8 in 
the presence of the applied voltage 

£ = {V\H |*) + <*| (1) 

subject to the following constraints: 

(ii) for the electrons flowing from electrodes into the device, the Wigner function at the 
two device-electrode contacts determined from \l/ in the presence of the voltage bias be equal 
to that obtained from the many-body ground state \l/o of the cluster in the absence of the 
bias; 

(iii) if not automatically satisfied,— 1 ^ the steady-state electric current J should be con- 
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strained to be position (q) independent 

(*|j ff |*) = J, (2) 

as required by the equation of continuity (dp/dt + V j = 0); 
(iv) the steady-state wave function is normalized 

(* l^> = L ( 3 ) 



In our recent work (Ref. l29f ). we showed that this approach is incorrect, as it fails to 
reproduce well-established results both for uncorrelated and correlated systems. The failure 
for the uncorrelated case is important for two reasons. First, in this case very large systems 
can be treated by this method without any other approximations than the method itself. 
This demonstrates that the completely unphysical results it generates are due to the method 
itself and cannot be assigned to the fact that, realistically, only very small parts of electrodes 
can be accounted for within accurate many-body calculations. Second, it challenges the need 
for and/or the usefulness of the Wigner function. Imposing boundary conditions by means 
of the Wigner function is only one ingredient of this method, and one may ask whether 
employing other boundary conditions could mender this variational approach. 

In the present paper, we shall consider such alternative choices. For the moment, we do 
not specify the properties to be constrained. Rather, instead of the above condition (ii), we 
shall consider general constraints of the form 

(V\Q K \V) = (y \Q K \y ) (4) 

for a set of Hermitian operators {Q K }- 30 Condition (ii) represents a particular case of Eq. 
flU), where Q K is chosen to be the Fano operator— 

F(x,p) = Y,aU a a x+l , a e~ 2 ^ h , (5) 

l,cr 

with a X)Cr and a\. a denoting the annihilation and creation operators for an electron of spin 
a located at site x. 

For the sake of simplicity, we shall consider one-dimensional discrete systems, with left 
(L) and right (R) electrodes containing electrons described within the tight-binding nearest- 
neighbor approximation (e = L,R). Although more general cases can also be considered, 
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we shall assume a total Hamiltonian of the form 
H = Hl + Hr + Hd + Ho, e 

H L = fX L a L a l,a ~ X K° { a la a l-l,a + a Ll,* a l,a) , 

Hr = Hr X a t,a a r,a ~ X V (<4,<r a r+l,a + 4+1,^^) , (6) 
r>qn,o- r>q R ,a 

^ D ,e = _ XX^ ( a l L ,<T a q L +l,a + a l L +l,a a q L ,o) ~ X] t lR^ {^R-h^qR,^ + a q R ,a a q R -l,a) 
a a 

-^£> = X/ X £ q a q,a a q,a ~ X X ^<?> CT (^.o^q+l.CT + a g+l,o- a g,o-) + ^int (^gt + l.er, • • • , n q R -l,a) ■ 
a <?=<?L + 1 °" <?=9L + 1 

Above, a; )0 .,a r)(T , and a gj0 - (0^,0^, and aJ jCr ) represent annihilation (creation) operators 
of electrons with spin a at sites /, r, and q belonging to the left (L) and right (R) 
electrodes, and to the device (D), respectively. fi LjR denote electrode chemical poten- 
tials, t's are hopping integrals, and e's are on-site energies in the device. The above 
model incorporates a rather general electron-electron interaction H int in the device. We 
need not to specify an explicit form of Hi nt , we solely assume that it can be expressed 
in terms of the electron numbers n q ^ = a qa a q ^ c . For concreteness, one can consider a 
particular explicit form H int = T l a l ,a 2 T iqL < qi , q2 < qR U qu a 1 ,q 2 ,a 2 n qiiai n q2j(T2 . Important par- 
ticular cases thereof are multi-site nanodevices described by extended Hubbard models 
,ai,q 2 ,a 2 ,q 2 $oi,-c2-> U qi! ai, q2 ,a 2 — ^<?i,<?2±i); or devices consisting of a single site 

(<?d), like point contacts or single-electron transistors (Anderson impurity model), which 
amounts to choose U qD)CriiqDiCr2 = U5 ai - a2) where U = or U 7^ 0, respectively. One should 
emphasize that the model specified above includes both uncorrelated (U qitaitq2tCT2 = 0) and 
strongly correlated cases. Phenomena like the Coulomb blockade or the Kondo effect de- 
scribed by the Anderson impurity model are typical examples of strong electron correlations, 
where the single-particle picture completely breaks down. 

Although not explicitly needed for the subsequent considerations, we give for illustration 
the form of the term W pertaining to the applied bias V 

V Qr1 V 

W= ^ X Vl,a + X X V A,° a q,° ~ ~o S a i,* a r,o, ( 7 ) 
l<q L ,a q=qL + l r>q R ,a 

where the concrete way how the potential V q drops across the device (— V/2 <V q < V/2) is 
not needed. 
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For model the electric current operator j q at site q (e = h = 1) has the expression^ 

cr 

If the total cluster possesses N sites, there are iV — 1 values (\I r |j g |\E r ). Therefore, the equation 
of continuity, Eq. 02]), yields N — 2 constraints of the form 

(*\j q -j qa \*) = 0, (9) 

where g stands for a fixed arbitrary site. For the moment, we need not to specify the Q K 's 
in Eqs. ©. 

To simplify the analysis, we shall assume that \l/o is nondegenerate. The wave function \I/ 
will be expanded in terms of the complete set of orthonormalized eigenstates of H (H\ty n ) = 

|*)= A |^ ) + E^n|^n). (10) 

n^O 

To simplify the discussion, we shall suppose that the model parameters entering Eq. (J6]) 
are real. This enables us to choose the eigenstates *ff n as real. However, the expansion 
coefficients A n can still be complex, and thus the general boundary conditions (jl]) and the 
equation of continuity (j2J) can be satisfied. We shall suppose that A Q is real, which amounts 
to choose the phase factor appropriately. 

Starting with a system in the ground state ^o, we shall look for the solution which 
minimizes the quantity 

(tf | H + W |*>-u;(tf |tf>-£ X K «*| Q K |*> - (* | |*o»+« E X* (*| l*> 

(11) 

with the supplementary constraints (EJ) , (HJ) , and (jHJ). This yields a set of equations, where- 
from the optimum values of the coefficients A and A n , and the (real) Lagrange multipliers 
u, Xqi ari d can t» e determined. 



III. LINEAR RESPONSE APPROXIMATION 



Similar to Ref. |29|, we shall work out the linear response approximation of the method 
described above, which should enable us to compute the zero-bias conductivity. Therefore, 
we shall only consider changes to the relevant quantities of the order 0(W). 



While in general this minimization represents a difficult nonlinear problem, it considerably 
simplifies in the linear response limit, applicable for a small applied bias. In this limit, 
the minimization amounts to solve a linear system of equations, which possesses a unique 
solution A = 1 + 0(W 2 ), A n = 0{W) for n ^ 0, \ K = 0(W), u = E + 0(W), and 
X q = 0(W). 

The quantities entering the minimization problem within the linear response approxima- 
tions are 

(V n \Q K \%) = Q K ,n = X K ,n + iy K ,n, (12) 

(tt n |W|*o)=W„, (13) 
(y n \j q -j qo \V )=iJ q ,n- (14) 

Above and throughout, the calligraphic symbols denote real matrix elements. The min- 
imization yields the following results in the linear response approximation. The expansion 
coefficients are expressed by 



"W n + E KQq,n + i E XqJq,n 

K q+qv 



The normalization condition fl3]) leads to u — E , and the constraints (Q and (ffl yield 



(15) 



£ n#0 (A n Q K , n + A n Ql n ) =0, (16) 

E n ^o (Ai — A*) Jq t n — 0. (17) 



At this point, it is useful to separate the real and imaginary parts of Eqs. ( ITojl . (jTBj) . and 
([17]) . This immediately leads to 

-V 2^ p p v l^Xql^^ ^r-Z^^ ^r^<z,n> 

Ex,' E P^f + E A.' E = 0. (19) 

Let us assume that there are Nl and Nr constraints of the form (jlj) imposed for the left and 
right electrodes, respectively. Eqs. ( IT8l) and ( TT9l) represent a linear set of + Nr + N — 2 
equations. Except for accidental cases where the determinants vanish, 32 these equations 
determine the N L + N R values X K and the N — 2 values x<? (<? <?o) °f the Lagrange 
multipliers uniquely. Once they are known, the expansion coefficients A n can be obtained 



from Eq. f[T5l) , which, in turn, allow to determine the q- independent steady state current 
J == Jq as 

J = ~ 2i E E A ^,„ + E x^n) (*« l^l *o) • (20) 

Because the angular parenthesis in the r.h.s. of the above equation is purely imaginary [see 
Eq- (E])], the quantity J is real. Notice that W« enters linearly Eqs. ( TT8l) and ( 1191) . and 
therefore the current computed from Eq. ( 1201) is proportional to the applied bias V. This 
means that the above minimization procedure yields the solution corresponding to the linear 
response limit. 

We end the part devoted to general considerations by addressing a technical issue. From 
the perspective of Ref. y, the variational approach discussed here would be ultimately 
intended to be used in conjunction with ab initio quantum chemical calculations for a real 
system, which comprises the device and parts of electrodes. It would be desirable that 
the latter are sufficiently large, such that the bulk electrode properties are approached. 
In practice, the size (N) of the system that can be investigated is inherently finite, and 
therefore an iV- dependence of the results is unavoidable [cf. Eq. (jSJ)]. The best one can hope 
in a realistic ab initio calculation is to be able to increase N until results converge, in a 
way similar to the much less demanding time-dependent DMRG calculations.— _ — In fact, 
increasing the size is so prohibitive that the most ambitious ab initio calculations can at 
most include a few electrode layers in the cluster used for transport calculations, and the 
saturation with increasing N cannot be systematically checked.- 33 Similarly, the expansion 
f fT0|) cannot exhaust the multi-electronic Hilbert space of a real system, which is infinitely 
dimensional. What one has to do there is to increase the number M of multi-electronic wave 
functions {^ n }n=\ until reaching convergence. Of course, the latter limitation only applies 
to a real system. For the discrete cases described by Eq. (jBJ), Af rapidly grows with iV and 
can become very large,— but it remains finite for any finite cluster. 

IV. CONSTRAINTS IMPOSED TO THE SINGLE-ELECTRON DISTRIBUTIONS 

24 . resorting to the Wigner 



So far, we did not specify the boundary conditions (j3j). In Ref. 
function (in a theoretical approach of the transport in correlated systems) was motivated 
by the fact that working with a correlated many-body wave function \I> precludes a descrip- 
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tion in terms of wave functions and energies e» for independent electrons. This has been 
interpreted as the impossibility of using single-electron Fermi-Dirac distributions /(cj) in a 
transport theory devoted to correlated systems.— This assertion might seem true, but it does 
not necessarily apply. As is well known, electron Fermi distribution functions are ubiqui- 
tously employed in transport theories, ranging from the (semi) classical Boltzmann equation 
to the Keldysh NEGF formalism. Whether for macroscopic, mesoscopic or nanoscopic sys- 
tems, whether applied to uncorrelated or strongly correlated systems, these theories have in 
common the assumptions that: 

(a) a separation in device and electrodes is possible, implying sufficiently small device- 
electrode couplings, which insignificantly perturb the electrodes. The properties of the 
electrodes connected to the device do not differ from those of the isolated electrodes; 

(b) electron correlations in electrodes are negligible, which implies that, there, the electron 
distributions are Fermi functions. 

Notice that only the electron distributions in electrodes are constrained, more precisely, 
they are constrained to be the same as in the isolated electrodes, and there they are Fermi 
functions irrespective whether the device is correlated or not. The Fermi distribution is the 
correct boundary condition for the uncorrrelated case. Because, as it will be shown below, 
with these boundary conditions the DG approach yields unphysical results, this suffices to 
demonstrate that the variational approach itself is incorrect. We do not intend to rigorously 
prove here that the Fermi distribution is the correct boundary condition for correlated 
nanoscopic/molecular systems. Still, we note that with boundary conditions expressed by 
Fermi functions one can successfully describe the electric transport in macroscopic systems. 
Moreover, by employing the same boundary conditions it was possible to explain delicate 
aspects of the transport in correlated nanosystems (e. g., the unitary limit for the Kondo 
plateau in single-electron transistors) in good agreement with experiment. Therefore, we 
argue that the description in terms of Fermi functions of the electrons in electrodes is 
plausible even if the latter are coupled to correlated devices. 

Let us now express the boundary conditions (HI) just in the aforementioned manner, 
i. e., by using the Fermi distribution. In our case, condition (a) requires sufficiently weak 
couplings t qL . and t qR)a between device and electrodes. In the absence of a voltage bias, 
the reservoirs are in equilibrium among themselves (/xj, = I^r = /•*)• An applied bias V 
does not affect the fact that each reservoir remains in local equilibrium, but it drives them 
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out of equilibrium with respect to one another, because it shifts the chemical potentials, 
A*l = ^ + V/2 and fiR = \i — V/2. It is just the imbalance of the corresponding chemical 
potentials that is kept constant by an external power supply, which causes a steady-state 
current through the device. 

Let us assume that for V = the single-electron states in the left electrode, which 
are specified by the labels k\,k,2,---, have the energies e Kl , e K2 , (Obviously, the con- 
siderations presented below also apply to the right electrode.) They can be obtained by 
diagonalizing the term Hl of Eq. 

H L = J2^ 3 ai 3 a Kj , (21) 
j 

although for the present purposes we need not to specify the explicit transformation T 

a la = ^T K]M a K . (22) 

j 

In the presence of a bias, the single-electron states remain specified by the same labels 
Ki, K2, ■ ■ ■■ it is the same linear orthogonal transformation T that diagonalizes Hl both for 
Hl = H and for n L = \i + V/2. The only change brought about by the bias is that the 
single-particle energies are simply shifted by V/2 with respect the former, Hl = Y^j{ e K,j + 
V/2)a\.a Kj . The occupancies of the single-particle states f Kj = (a\.a Kj ) at equilibrium 
{V = 0) and in the steady state {V ^ 0) are the same: f^ R (V) = f^ R = f^ R (V = 0). 
Only their energy distributions are shifted f L ' R (s K .) = fo ,R {s Kj =F V/2). The situation is 
schematically depicted in Fig. [TJ Of course, the above considerations are not restricted to 
the linear response limit. 

By choosing now the electron occupancies to express the boundary conditions (Q K = 
aj.au,) and assuming, as usual, that the electron distributions in the reservoirs are not affected 
by the coupling to the device, we can reexpress Eqs. (jl]) as 

(tf|aia«|tf) = <tt |4a*l*o), (23) 

which apply both to the left and the right electrodes. Most importantly, the quantities 
Qn,n — Xn,n = (^n\ ot\.a K |\I/o) entering the l.h.s. of Eq. (TT2"j) are real. With this choice, 
y K}n = 0, and Eqs. ( fTBi) and ( |T9l) become two sets of decoupled equations for A's and x' s > 
respectively 

A «' 1^ W~ = 2^ p ET-V, 1/4) 

10 




FIG. 1: (Color online) Effect of an applied bias V on the single-electron distributions / and f R in 
the left (L) and right (R) electrodes. At V = 0, the electrodes are in equilibrium, and their chemical 
potentials are equal, e r f = e r f = (black dashed horizontal line). The bias V creates an imbalance 
between the chemical potentials fJ-L = £ k f + V/2 and Hr = £ k f — V/2 responsible for the current 
flow. The figure visualizes the fact that the bias only shifts the single-electron energies e KL R in the 
two electrodes by ±V/2 in opposite directions: e n % (V) = e R % + V/2 and e K i (V) = e K i — V/2. 

L L ,R Ft Ft 

However, the bias affects neither the labels of the single-electron states [k 1 l r(V) = k 1 l r ] nor their 
occupancies: f%* (V) = f L > R (e K i ± V/2; V) = f L > R (e K i ; V = 0) = f^ R (V = 0) [cf. Eq. ©]. 

K L,R L,R L,R L,R 

The representation in the figure is schematic: in reality, neither the energies e (vertical axes) nor 
the quantum indices kl,r (oblique axes) have equidistant values. 
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Ex 9 'E^f = 0. (25) 
From the latter equations we get \q — 0,— and then Eq. (120]) yields 

J = 0. (26) 

This means that the approach considered above predicts a vanishing zero-bias conductance 
(G = J/V = 0) at least for the rather general models of the type One should note 
that this prediction also holds if not all k's, but only some of the occupancies of the single- 
electron states are subject to constraints of the type (j23|) . Furthermore, the prediction (126]) 
also applies to the case where the constraints ([9]) are not applied for all the sites q, but, e. g., 
only at the two contacts (i. e., the current flowing from the left electrode into the device is 
equal to that from the device into the right electrode). 

Needless to say, the prediction fl26l) of a vanishing current in the linear response limit 
(vanishing zero-bias conductance) is completely unphysical and flagrantly contradicts nu- 
merous well established results. The existence of a nonvanishing zero-bias conductance 
G = Gq = 2e 2 /h in the textbook example of conduction through a resonant single level 
system,- or the fact that the same value Go characterizes the Kondo plateau in the transport 
through a single-electron transistor— 1 ^ are only two particular cases, which are described 
by the general model ([6]) for which the above considerations apply. 

Concerning the linear response limit of the variational approach discussed here, we still 
have the following comments. Without imposing any boundary conditions of the type (jlj), 
the minimization of £ of Eq. (CQ) yields equations identical to Eqs. (fT5l) and ( JT9l) if we 
set A K = 0. Because in this case the latter group of equations only possess the trivial 
solution Xq = 0, the expansion coefficients obtained from Eq. ( fl5l) have the form A n = 
A n = W n /(E — E n ), i. e., the solution |^) of Eq. fflOl) is nothing but the ground state \Q) 
of the total Hamiltonian H + W obtained within the first-order perturbation theory with 
respect to W . Being the ground state, that is, an eigenstate of the total Hamiltonian, it 
obviously obeys the equation of continuity: the current is site independent, more precisely 
J = J q = (Q\j q \Q) = 0. 3 - When the supplementary (boundary) conditions (HI) are imposed, 
the corresponding Lagrange multipliers become nonvanishing X K ^ 0, as seen from Eq. (f!8l) 
or even from the more particular Eq. (f24l) . Consequently, A n ^ A n , that is, the solution of 
optimization does differ from the ground state with applied bias, |^) ^ \Q). The essential 
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point is that, as we have seen above, in order to sustain a nonvanishing current (J ^ 0), the 
matrix elements of the operators employed to impose the boundary conditions should have 
nonvanishing imaginary parts, lmQ Kn = y K>n 7^ 0. Luckily, this happens to be the case if the 
Fano operator §5§ is used in Eq. (J3J), which amounts to formulate the boundary conditions in 
terms of the Wigner function.— This choice yields nonvanishing currents, which sometimes, 
by chance, could be comparable with measured values^ 1 ^ or look plausible.— 1 ^ However, 
as we have recently unambiguously demonstrated,— the variational approach based on the 
Wigner function fails to recover the simplest results both for uncorrelated and correlated 
nanosystems, where its predictions are completely unphysical. 



V. DISCUSSION AND CONCLUSION 



In Ref. 123, we demonstrated that the DG approach^— fails to reproduce the simplest 
well established results for transport in nanosystems and especially criticized the boundary 
conditions. In view of the fact that a reliable method to deal with correlated electron 
transport is hardly needed, one might still think to be able to develop the valid approach 
based on the DG method but employing other boundary conditions. Therefore, in the 
present paper, we have inquired the validity of the DG variational approach by considering 
constraints other than those imposed to the Wigner function at the boundaries. We believe 
that this investigation is important for specialists in the field thinking to apply such a 
modified DG scheme: the implementation for realistic ab initio calculations is by no means 
an easy task, and this probably explains why no other group applied the DG method in 
spite of its claimed success. 

In the first part we have examined the linear response approximation by considering 
general boundary conditions. Further, we have critically analyzed the claim of DG on 
the Fermi-Dirac single-electron distributions. Namely, they claimed^— that electron Fermi 
distribution functions, which are ubiquitously employed in transport theories, cannot be used 
in correlated systems because the single-particle picture breaks down, and resorting to the 
Wigner function was proposed as a way out of this difficulty In Sect. IIVI we have explained 
that strong electron correlations in devices do not preclude the usage of single-electron 
distributions to express boundary conditions. The reason is that they should be imposed 
in electrodes, and there electrons are (assumed to be) uncorrelated. In fact, constraining 
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single-electron distributions to be of Fermi-Dirac type in electrodes is the common feature 
of the most transport theories employed for macroscopic, mesoscopic and nanoscopic scale, 
and of interest there is just the case of correlated systems. 

Obvioulsy, it is possible to impose constraints on the single-electron distributions in un- 
correlated systems. As an important point of our analysis, we have demonstrated that, in 
particular, even for such uncorrelated systems and these standard boundary conditions (a 
case beyond of all question) the variational approach predicts a vanishing zero-bias conduc- 
tance, which represents a completely unphysical result. Although this very fact is sufficient 
to invalidate this variational approach, the conclusion of the present investigation applies 
to more much general situations where strong electron correlations are present. Such more 
general examples include (but are not limited to) the systems described by the model Hamil- 
tonian (jfjj). 

In principle, any theoretical approach of transport is based on two main ingredients, (i) 
First, one has to determine a certain quantity, which characterizes the transport cluster, (ii) 
Second, one has to impose boundary conditions to link the cluster to electrodes. The first 
ingredient is e. g., the transmission coefficient, Green function(s), or transition probabilities 
in Ladauer, NEGF, or master equation approaches, respectively. While this first ingredient 
is different from one transport approach to another, the second is common for virtually 
all approaches: whether correlated or not, the cluster is linked to electrodes by employing 
Fermi distributions. Because none of these approaches yields results which are manifestly 
unphysical, one has to admit that the formulation of the boundary conditions in terms of 
Fermi distribution functions is legitimate, or at least does not lead to manifest absurdities. 
If these boundary conditions were wrong, all these widely employed approaches would also 
yield unphysical results, but this is not the case. In fact, it is hard to imagine a physical 
situation more typical for using single-particle Fermi distributions than to describe electrons 
in metals (electrodes). 

The approach of Delaney and Greer also comprises these two aspects: (i) One minimizes 
the total energy E of the finite cluster used for transport calculations in the presence of 
an applied voltage. The energy of this open system is computed as the average value of a 
Hermitian Hamiltonian, E = (H + W). (ii) This is a constrained minimization: specifically, 
boundary conditions are applied to the Wigner function. 

In Ref. |29|, we demonstrated that, in the form proposed in Ref. |24|, this approach yields 
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unphysical results and especially criticized the boundary conditions adopted by Delaney and 
Greer in the context of their approach. In the present paper, we have lifted these specific 
constraints and replaced them with the widely used boundary conditions, expressed in terms 
of Fermi functions. The result obtained, a vanishingly zero-bias conductance, is again quite 
unphysical. So, even with these usual boundary conditions, formulated in terms of Fermi 
distributions in electrodes, in the way common to the most widely employed transport 
theories, the approach based on the DG- variational method fails to pass the minimal decisive 
test, which any approach of electric transport at nanoscale must satisfy; namely, to be at 
least able to correctly describe the conductance in the linear response (Kubo) approximation 
or even to reproduce well established results of the Landauer theory in the absence of 
correlations (see, for example, Ref. l38l ). 

The main result of our study is that the DG- variational approach does not work even 
with modified boundary conditions. This implies that the DG-approach is invalid not (or not 
only) because of the boundary conditions. If the imposition of the boundary conditions in 
terms of Wigner function were the only wrong point, the variational approach would produce 
correct results with correct boundary conditions. Doubtless, the Fermi functions represent 
the correct boundaries for uncorrelated system, a fact which proves the failure of the DG 
approach for that case. If the DG variational approach is incorrect even for uncorrelated 
systems, it is inconceivable that it works for correlated systems. We are then necessarily 
led to the conclusion that the other ingredient (i) of the DG- variational approach, namely 
the minimization of the cluster's total energy, computed as the average of a Hermitian 
Hamiltonian, is inappropriate. 

So, the failure of the approach of Delaney and Greer cannot be solely assigned to the 
manner in which the boundary conditions are imposed. Constraining the Wigner function 
f{q,p) at the boundaries can also produce valuable physical results for the steady state 
current, provided that f(q,p) is determined from the Liouville equation for df/dt = and 
the values f(qL,R,p) at the boundaries determined by the Fermi functions in reservoirs.— 
However, in that case, the fact that the investigated systems are open leads to a non- 
Hermitian Liouville operator, which yields time irreversibility and currents that saturate to 
values characterizing the steady state flow. The counterpart of this procedure would be to 
employ non-Hermitian Hamilton operators H 7^ %\ amounting to consider imaginary self- 
energies resulting from electrode-device interactions.-- Attempting to develop a variational 
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approach for a finite open system (device) described by a non-Hermitian Hamiltonian, e. g., 
by minimizing {CM. + W — E) + W — Eyj instead of Eq. (pQ), might be an interesting 
alternative to the existing approaches to correlated transport, but to our knowledge such 
a method is not available at present. Till then, one is forced to admit that the variational 
procedure attempting to obtain the steady state current through a finite open system de- 
scribed by a Hermitian Hamiltonian and a wave function that minimizes the total energy in 
the presence of a voltage bias with certain boundary conditions is unable to reproduce well 
established results, on which there is an incontestable agreement between experiments and 
other theoretical approaches. 
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